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Microlensing of an extended source by a power-law mass 
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ABSTRACT 

Microlensing promises to be a powerful tool for studying distant galaxies and quasars. 
As the data and models improve, there are systematic effects that need to be explored. 
Quasar continuum and broad-line regions may respond differently to microlensing due 
to their different sizes; to understand this effect, we study microlensing of finite sources 
by a mass function of stars. We find that microlensing is insensitive to the slope of the 
mass function but does depend on the mass range. For negative parity images, diluting 
the stellar population with dark matter increases the magnification dispersion for small 
sources and decreases it for large sources. This implies that the quasar continuum 
and broad-line regions may experience very different microlensing in negative-parity 
lensed images. We confirm earlier conclusions that the surface brightness profile and 
geometry of the source have little effect on microlensing. Finally, we consider non- 
circular sources. We show that elliptical sources that are aligned with the direction of 
shear have larger magnification dispersions than sources with perpendicular alignment, 
an effect that becomes more prominent as the ellipticity increases. Elongated sources 
can lead to more rapid variability than circular sources, which raises the prospect of 
using microlensing to probe source shape. 
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1 INTRODUCTION 

Microlensing is an increasingly important tool for study- 
ing small-scale structure in lens galaxies and source quasars 
In recent years, microlensing has bee n observed in a num- 
ber of multiply - lensed quasars (e.g.. IWozniak et alj l200d: 



Schechter et al 



20031: IRichards et al.1 12004 IKeeton et al l 
20061 : iParaficz et al.l 120061 ) . Mi crolensing modeli ng has also 
been improving. For example. iKochan ek (2004) has intro- 
duced a sophisticated technique for analyzing light curves. 
Even so, there are aspects of the models for the lens and 
source that still need to be considered for microlensing to 
reach its full potential. This is especially important for 
quasar microlensing, where several length scales are in- 
volved. 

Different regions of a quasar emit radiation in roughly 
distinct bands. For example, continuum (blackbody) radia- 
tion in the optical and x-ray bands is emitted from the ac- 
cretion disk surrounding the supermassive black hole, while 
broad emission lines in the optical and UV are thought 
to originate from clouds in a region outside of and larger 
than the accretion disk. Radio emission comes from even 
larger structures. Roughly speaking, a source can only be 



affected by objects in the lens galaxy whose Einstein radii 
are comparable to or larger than the source size. Conse- 
quently, the continuum, broad-line, and radio regions probe 
different structures in the lens galaxy. Radio jets can be 
used to probe dark matter su bstructure (jMetcalf fc MadatJ 
l200ll : [Dalai fc Kocha nck 2002 ) in lens galaxies while the ac- 
cretion disk (jjaroszvnski. Wambsganss fc Paczvriskil Il992l: 
Mortonson. Schechter fc Wambsganssl 120051 ; iPoolev et al.l 
120061 ) and broad-line regi on (BLR) are used for study- 
ing the stellar component (ISchneider fc Wambsganssl [l990l : 
IRichards et alj|2004l : IKeeton et alj|2006l ). In principle, both 
methods can also be used to examine the light source. This 
is of particular interest fo r the BLR whose prope rties are 
not well understood (e.g., IPeterson fc Horndliooa ). In this 
paper we investigate the potential of microlensing to deepen 
our knowledge of the BLR and accretion disk and to deter- 
mine both the abundance and mass function of stars in lens 
galaxies. 

Until recently, sources relevant for microlensing could 
not be resolved. Therefore, many theoretical models have 
assumed a point source, and have focused on deter- 
mining properties of the lensing galaxy. Such investiga- 
tions have found that microlensing magnification distri- 
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butions are not very sensitive to the shape of the mi- 
crolens mass function if i t spans an order of magni- 
tude or so in mass (e.g., IWvithe fc Turner! l200lh . The 
magnification distributions do look different if there are 
two distinct mass components: not just stars, but also 
"dark matter" that coul d come in the form of a s moot h 
mass component (e.g., ISchechter fc Wambsganssl [2002), 
or a set of discre te objects that are much less mas- 
sive than the stars (Schcchtcr, Wambsgan ss fc Lewis! [20o3 ; 



iLewis fc Gil-Merindl2006h . We generalize the previous stud- 
ies by considering microlensing of an extended source. The 
source size Rs introduces a new length scale, which heuristi- 
cally divides the microlenses into two categories: microlenses 
with Re > Rs are massive enough to be felt individually, 
while microlenses with Re < -Rs act like a smooth com- 
ponent. We therefore conjecture that the finite source size 
makes the magnification distributions sensitive to the mi- 
crolens mass f unction even when it spans just 1-1.5 orders 
of magnitude. ILewis fc Gil-Merinol (2006) recently studied 
microlensing of an extended source by a bimodal mass func- 
tion; we now consider a continuous mass function. 

Studying microlensing of an extended source is espe- 
cially relevant for the BLR, because recent o bservations sug- 
gest i t has a size Rblr ~ 10 16 10 18 cm (I Richards et al.l 
120041 : iKeeton et al.l l2006t ) comparable to stellar Einstein 
radii. Understanding the effects of source size should help us 
probe BLR length scales more precisely. We perhaps should 
not expect to probe the BLR surface brightness di s tribu - 
tion, however: iMortonson. Schechter fc Wambsganssl (|2005l ) 
suggest that microlensing is not very sensitive to the source 
brightness profile. Their assertion is based on simulations of 
circular sources with several specific surface brightness pro- 
files. We perform similar "numerical experiments" to con- 
sider other source properties, notably shape and geometry. 

The possibility of a non-circular source has not been 
considered in previous microlensing analyses. However, ac- 
cretion disks viewed at random inclinations would generi- 
cally appear elliptical rather than circular. A similar situa- 
tion may apply to the BLR if i t has a disky structure (e.g ., 
iMurrav fc Chiang! 1 19981 : lElvisI |2000| ; iRichards et~aH 12004) . 
We also consider annular accretion disks. Such models are 
important for two reasons. First, quasar accretion disks have 
inner radii defined by the innermost stable circular orbit of 
a particle in motion around the central black hole. Second, 
typical models (e.g., the Shakura-Sunyaev disk), emit over 
a wide range of wavelengths, with each band corresponding 
to different annular regions. 

This paper is organized as follows. Lens modeling and 
simulations are discussed in jj2] Our results are presented in 
Sj3] We first consider the general problem of how source size 
and lens mass impact microlensing (333}- In following sub- 
sections, we investigate the effects of dark matter ( £I3.2[) . 
source profile ( i]3.3|l . source ellipticity and position angle 
(" i]3.4|l . and accretion disk geometry (Sj33}. We construct 
light curves in i )3.6l to investigate variability timescales. Our 
conclusions are summarized in 34] 



in the lens galaxy around a lensed image. The size of the 
region is chosen to satisfy two conditions. First, it must be 
large compared to a typical stellar Einstein radius, which is 
the relevant scale for microlensing. Second, it must be small 
compared to the scale of the lens galaxy, namely the image 
separation. The latter allows us to take the mean densities of 
stars and dark matter to be constant. These criteria are not 
very restrictive since Einstein radii are typically on scales 
of microarcseconds while image separations are on scales of 
arcseconds. 

We describe the stellar population of the lens by a mass 
function dN/dm, which gives the number of stars per unit 
mass between m and m + dm. We use a power law of the 
form 



dN 
dm 



(mi ^ m ^ 7712) 



(1) 



for some mi and m.2. Rather than specifying the mass limits 
mi and mi explicitly, we adopt the equivalent approach of 
giving the ratio m 1 /m2 along with the mean mass m. A typ- 
ical choice for the power law index is a — 2.35, the Salpeter 
initial mass function. The mass function is normalized so 
that the scaled mass density is re, . In addition to stars, the 
galaxy may include a continuous component with density 
re c . The final ingredient is the shear 7, which accounts for 
tidal forces from outside the patch of stars being considered. 
To obtain a relation between re and 7 we must choose a mass 
model for the lens galaxy. We use a singular isothermal el- 
lipsoid for which k = 7. This model is si mple and provides a 
reasonable fit to observed systems (e.g..lTreu fc Koopmand 
I2004 I Rusin fc Kochanekll2005l : iTreu et alj|2006h . 

In the absence of microlensing, the magnification is 
given by 



Ho = [(1 - re) 2 - -y 2 Y 



(2) 



where We consider a typical bright image with 

(io = ±10 corresponding to re = 7 = 0.5=f0.05. In microlens- 
ing, the spatial distribution of stars is random, so the mag- 
nification at a given time will be drawn from a probability 
distribution. If the distribution is broad, then chances are 
high that the magnification will differ from that predicted 
for a smooth model, viz. ^o- The effects of microlensing are 
therefore naturally described by the dispersion (standard 
deviation) of the probability distribution. 

Computing the magnification distribution for given re, 
7 and dN/dm can o nly be done n umeri cally. We use the 
microlensing code of IWambsganssI (|l999h , which gives the 
magnification of a point sourc43 as a function of position 
over a square region of the source plane with side length 
L — 15Re{Mq). In mean mass units the side length is given 
by L = 15i?s(m) for mi/m^ = 1, L = 31.71i?E(m) for 
mi/7712 = 0.1, and L — 52.26i?s(m) for mi/ma = 0.03. We 
create magnification maps with a resolution of L/1024 (see 
Figure [TJ . To obtain a statistical sample, we perform 100 
realizations for each set of parameters we consider. 

We must convolve the magnification map with a surface 
brightness profile in order to compute the magnification of 



2 METHODS 

We consider microlensing of an extended source by a distri- 
bution of stars and dark matter. We zoom in on a region 



1 Strictly speaking, the map gives the magnification of a source 
with the shape of the pixel. In practice, the source sizes that 
interest us are much larger than the pixels. 
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Figure 1. Magnification maps for a positive-parity image with k = 7 = 0.45, implying a macro-magnification of /io = fO. The gray 
scale indicates the magnification with values in the range fi = f (black) and fj, = 30 (white). Panels show Salpeter mass functions with 
mi/m2 = f (left), O.f (middle), and 0.03 (right). Each map has a side length L = 15iiE(m). Magnification histograms are generated by 
convolving the surface brightness profile of the source (indicated by circles) with the magnification map. 



an extended source. We use Gaussian, uniform and linear 
profiles, respectively denned by 

R 2 In 2" 



for radii satisfying 



h(R) 
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(3) 



(4) 



(5) 



where Rs is the half-light radius of the source, and the 
sources are normalized to unit flux. In microlensing, the nat- 
ural l ength scale is the Ein stein radius of the mean mass star 
(e.g., iLewis fc Irwinlfl996l ). We henceforth quote the source 
size as rs = Rs /-Rfi(m). 

These models are simple but nevertheless use- 
ful. The Gaussian model is popular in microlens- 
ing s tudies (e.g., Wambganss. Paczvhski fc Schneider! 



199Q|; IWvithe. Agol fc Fluke! l2002h . so it is good to 
include. The uniform disk is the simplest model 
conceivable, while the linear disk has at least one 
physical c onnection. These same three mod els wer e 
used by iMortonson. Schechter fc Wambsganssl (|2005T ) . 
so we can compare their results with ours. While 
IMortonson. Schechter fcTW ambsg anssI (|2005l ) (and many 
others) applied the models to accretion disks, we imagine 
that they are useful for describing BLRs as well. In par- 
ticular^the linear disk is similar to the biconical BLR of 



lAbaias et all (|2002T ). 

We consider an elliptical source by making the substi- 
tution 



I(R) - I(p)/q, 



(6) 



with minor-to- major axis ratio q and elliptical radius defined 
by p 2 = R 2 cos 2 8 + R 2 sin 2 9/q 2 . We also consider a uniform 
annular disk, with inner-to-outer radius ratio Q, by making 
the replacement 



h(R) h{R) 



l + < 



(7) 



V2QR& 
l + Q 2 



V2R S 
l + Q' 



(8) 



In the following section our fiducial model assumes a 
Gaussian source and a lens described by a stellar popula- 
tion whose masses are given by a Salpeter mass function 
with m\/m2 = 0.1. We assume that k = ft, = 0.5 =F 0.05 
for images of positive and negative parity, respectively. We 
explicitly state when other models are to be used. 



3 RESULTS 

We now study microlensing of an extended source by a 
power-law mass distribution. IWvithe fc Turner! \ 2001) con- 
clude that point-source microlensing magnification distribu- 
tions are not significantly affected by the choice of mass 
function. We consider whether this result can be extended 
to the case of a finite source ( £|3.1[) . Dark matter can affect 
microlensing in surprising ways, raising the possibility of us- 
ing microlensing to measure the density of dark matter a t 
the image positions (|Schechter fc W ambsganss 2 0021 . [20041 ') . 
We generalize earlier work by including a mass function of 
stars and an extended source ( £|3.2ll . 

We also explore how varying properties of the source 
impacts microlensing magnification distributions. In H3.3I 
we describe a source by three surface brightness profiles. We 
broaden the discussion in £13.41 to include elliptical sources. 
Finally, we consider annular sources in 



3.1 Source Size and Lens Mass 

We begin by examining how microlensing of a finite source is 
affected when we vary t he mass range and l ogarithmic slope 
of the mass function (cf. IWambsgansslll992h . Figure [2] shows 
the mass functions we use. Figures [3] and [4] show magnifica- 
tion histograms for the different mass functions and different 
source sizes, for positive and negative parity. First consider 
the effects of the mass range, shown in the top row of each 
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Figure 3. Magnification histograms for different source sizes and mass functions for a positive-parity image with \fio\ = 10. Columns 
show source sizes in mean-mass Einstein radii of rg =0.1, 0.4, 0.7 and 1.0. The top row shows mass functions with logarithmic slopes 
of a = 2.35. The solid, dotted, and dashed curves have mi/rri2=l, 0.1 and 0.03, respectively (see Figure l2l top panel). The bottom row 
shows mass functions with mi/m2=0.1. The solid, dotted, and dashed curves have a =1.85, 2.35 and 2.85, respectively (see Figure [2] 
bottom panel). 
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Figure 2. Mass functions with fixed means, rh, for various mass 
ranges, mi/m2, and logarithmic slopes, a. The top panel shows 
mass functions with a = 2.35. The solid, dotted and dashed 
curves have mi/rri2= 1, 0.1, and 0.03, respectively. The bottom 
panel shows mass functions with mi/r?i2 = 0.1. The solid, dotted 
and dashed curves show a= 1.85, 2.35 and 2.85, respectively. 



figure. Increasing the mass range causes the magnification 
distribution to broaden slightly, especially for larger sources. 
This is shown more directly in the top panel of Figure [S] 
which plots the magnification dispersion versus source size 
for the different mass ranges (for the positive parity case). 
When the source is small we recover the previous result that 
the mass range does not affect the magnification distribution 
l|Wvithe fc Turner1l200ll ). However, as the source gets larger 
there is more difference between the three mass functions. 

To understand why the magnification dispersion in- 
creases as the mass range increases, we return to Figure [2] 
The top panel shows that increasing the mass range causes 
the mass function to "spread out" : the lower limit decreases 
slightly, while the upper limit can increase substantially. A 
high upper limit allows massive stars to exist, although they 
will be fairly rare because the mass function is steep. Thus, 
some magnification maps will contain one or a few massive 
stars that significantly affect the microlensing, while many 
will not. We believe this explains why increasing the mass 
range increases the magnification dispersion. It also explains 
why the mass range becomes more important as the source 
size increases: large sources are most sensitive to massive 
stars. 

Now consider the slope of the mass function. Figures 
[3] [5] show that the slope hardly affects microlensing at all, 
regardless of the source size. The reason is that changing 
the slope shifts the mass function left or right (see Figure 
[2]), but not dramatically. 
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Figure 4. Similar to Figure [3] but for a negative parity image with \/io \ = 10. 
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Figure 5. Dispersion in log \fi\ computed from histograms for 
source sizes in the range ^ rg ^ 4 (see, e.g., Fig. [3j versus source 
size for different mass functions, for a positive-parity image with 
\fio\ = 10. The top panel shows mass functions with logarithmic 
slopes of a = 2.35. The solid, dotted, and dashed curves have 
mi/m.2 = 1, 0.1 and 0.03, respectively. The bottom panel shows 
mass functions with mi/m2 = 0.1. The solid, dotted, and dashed 
curves have a = 1.85, 2.35 and 2.85, respectively. 



We conclude that microlensing of an extended source 
may offer the possibility of determining the dynamic range 
(7711/7712) of the stellar mass function, but not for deter- 
mining the mass function slope. In light of this result, we 
henceforth restrict our attention to a Salpeter mass function 
(a = 2.35), and we focus attention on the case 7711/777,2 = 0.1. 



3.2 Dark Matter Content 

We now consider how the mass fraction in dark matter af- 
fects microlensing. Controversy remains as to whether cos- 
mological simulatio ns of dark matter agree with galaxy 



;e with galaxy 
Havnesl 12005; 



observations (e.g., ISpekkens. Giovanelli 
iGerhardl l2006j and references therein). Part of the prob- 
lem is that most observations (galaxy dynamics, grav- 
itational macrolensing) depend on the global mass dis- 
tribu tion in a galaxy, rather than the l ocal mass den- 
sity. ISchechter fc Wambseanssl (|2002l . l2004h argue that mi- 
crolensing offers a new way to measure the density of dark 
matter at the positions of the lensed images. They consider 
a uniform mass function and a point source; we generalize 
the analysis to a mass function and an extended source. 

Figure [6] shows magnification histograms for dark mat- 
ter mass fractions of f c = n c /n = 0.2, 0.4, 0.6, 0.8, and 
0.99. One might expect that cLS Kic IS increased, microlens- 
ing would become less important, since fewer stars produce 
simpler caustic networks. For k c — > k the magnification dis- 
tribution indeed approaches a 5-function centered at fio, as 
seen in the bottom row of Figure [6] 

For smaller values of k c , however, the histograms show 
more structure. In particular, sec ondary peaks app ear in 
several of the histograms (see, e.g.. lRauch et al.ll 19921 ). The 
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Figure 6. Magnification histograms for an image with j/io = 10 
for different dark matter mass fractions. From top to bottom, 
panels show / c = 0.2, 0.4, 0.6, 0.8, and 0.99. The left and right 
columns have rg = 0.1, and 0.5, respectively. The solid and dotted 
curves have positive and negative parity, respectively. We use a 
Salpeter mass function (a = 2.35) with mxjm^ = 0.1. 



different peaks correspond to different numbers o f microim- 
ages (see lGranot. Schechter fc Wam bsganss 2003, especially 
Figure 4). If dark matter is the primary mass component, 
the probability that a source will have multiple microimages 
is low. In that case, the magnification distribution has a sin- 
gle peak near the expected magnification in the absence of 
microlensing (e.g., f c — 0.99). For smaller values of / c , the 
density of caustics increases, which in turn raises the prob- 
ability that a source will have extra microimage pairs. The 
magnification distribution therefore acquires a second peak 
associated with regions of the source plane for which an ex- 
tra image pair is created (see left-hand column of Figure [6}. 
When f c is low, regions with extra microimages become the 
norm, and the peaks corresponding to different numbers of 
microimages become less distinct. Also, an extended source 
often covers regions with different numbers of microimages, 
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Figure 7. Dispersion in log \fj,\ computed from histograms for 
source sizes in the range ^ rg ^ 5 (see, e.g., Fig. |6j versus 
source size for different dark matter fractions. The top (bottom) 
panel shows an image with positive (negative) parity. The solid, 
dotted, dashed, dot-dashed and dot-dot-dashed curves have f c = 
0.2, 0.4, 0.6, 0.8, and 0.99, respectively. 



smearing out the effects of additional microimages (see right- 
hand column of Figure [6} . 

Figure[6]also demonstrates the importance of parity. We 
find that negative parity i mages lead to distributions with 
tails at low magnification. ISchechter fc Wambsganssl (|2002D 
find the same behaviour for a point source. As the source size 
is increased, the tails become less apparent. When rg = 1.0 
(not shown), the two parities give nearly identical results. 
One surprise is that a difference between positive and nega- 
tive parity can be observed in the skewness of the distribu- 
tions even for f c — 0.99. This means that even a small stel- 
lar mass fraction gives rise to noticeable parity-dependent 
effects. 

Figure [7] uses the magnification dispersion to quantify 
the effects of parity and source size. In the positive parity 
case (top panel), replacing stars with dark matter decreases 
the dispersion for all source sizes, which makes intuitive 
sense. However, in the negative parity case (bottom panel) 
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when the source is small, diluting the stars with dark mat- 
ter increases the dispersion, at least in the range f c ^ 0.8. 
ISchechter fc Wambsganssl (|2002l l first found this result for a 
point source and a uniform stellar mass function. We now see 
that it holds for small extended sources as well. We discover 
though, that when the source is large, the trend reverses: 
increasing the dark matter fraction decreases the magnifica- 
tion dispersion. It seems notable that the curves of disper- 
sion versus source size for different dark matter fractions all 
cross at roughly the same source size (rs ~ 0.8), although 
we do not know whether this is significant. 

It is worth p ointing out that 

iDobler. Keeton fc Wambsganssl HoOff) examine the ef- 
fects of dark matter and source size on the magnification 
for demagnified lensed images. They find that increasing 
the dark matter fraction always decreases the magnification 
dispersion for a demagnified negative parity image. (Recall 
that our negative parity image is magnified.) However, for 
a demagnified central image, diluting the stars with dark 
matter increases the dispersion for a small source, but de- 
creases the dispersion for a large source. Direct comparison 
between those results and ours is not possible due to the 
different macro parameters. Nevertheless, it seems clear 
that the effects of dark matter on microlensing depend in 
a complicated way on the parity, the macro-magnification, 
and the source size. 

Our findings imply that the continuum and broad- 
line regions may experience very different microlensing in 
negative-parity lensed images. The continuum emission re- 
gion is small and should therefore have a magnification dis- 
persion that increases with the dark matter fraction. By 
contrast, in many cases, the BLR may be large enough that 
the disp ersion will decrease as the dark matter f raction in- 
creases (|Richards et al.ll2004l ; iKeeton et~ai] 12006). For posi- 
tive parity images, the continuum and BLR will both have 
a dispersion that decreases with the dark matter fraction. 
This may turn out to be a very important physical effect 
allowing us to probe both the dark matter content of lens 
galaxies and the structure of lensed quasars. 

3.3 Source Profile 

In the remaining subsections, we return to models consist- 
ing of a purely stellar mass component, and consider how 
microlensing depends on properties of the source. We first 
exami ne different source surface brightnes s profi les, follow- 
ing iMortonson. Schechter fc Wambsganssl (|2005h . Figure [8] 
shows the dispersion versus source size for Gaussian, uni- 
form, and linear profiles (defined in ^2]). We see that the dis- 
persion decreases as the profile becomes steeper, although 
the effect is not strong. We therefore confirm that the dis- 
persion depends weakly on the source profile. 

3.4 Ellipticity and Position Angle 

We now allow the source to be non-circular. This pos- 
sibility has not been considered in previous microlensing 
analyses, although it is an important physical effect (see 
iKochanek et al.ll2006r i. For a population of thin disks with 
random inclinations, a face-on source is rare; the average 
projected ellipticity is e = 0.5. Therefore, models of mi- 
crolensing need to allow a non-circular shape for the source. 



T 




1 2 3 4 5 

Source Size r s 

Figure 8. Dispersion in log versus source size for a posi- 
tive parity image with different source profiles. Solid, dotted and 
dashed curves have uniform, linear and Gaussian source profiles, 
respectively. A Salpeter mass function with mi/m.2 = 0.1 is used. 



This is certainly true for continuum microlensing, whose 
source is presumably a thin accretion disk. It may be true for 
broad-line microlensing as well, given e vidence that BLRs 
may have some disky structure (e. g., iMurrav fc Chiangl 
ll998l : lElvisll2000l : lRichards et al.ll2004h . 

Figure [9] shows how the magnification dispersion de- 
pends on the ellipticity, e = 1 — q, and position angle, PA, of 
the source. In each panel we see that the dispersion increases 
monotonically with position angle, an effect which becomes 
more pronounced for large ellipticities. To understand this 
behaviour, first note that PA=0° describes a source whose 
major axis is orthogonal to the direction of shear, which de- 
fines the long axis of the caustics. An extended source with 
PA=0° is likely to cover one or more caustics regardless of 
where it is centered (see Figure [TO)) . Small changes in the 
source position do not produce dramatic changes in the mag- 
nification. By contrast, for a source with PA=90° (aligned 
with the caustics), small displacements of the source can 
change the number of caustics covered, with corresponding 
large deviations in the magnification. This explains why the 
magnification dispersion is higher for sources aligned with 
the caustics (PA=90°) than for orthogonal sources (PA=0°). 
These effects become even more pronounced for more highly 
elongated sources. 

In Figure [11] most of the effects discussed so far are 
considered simultaneously. As in Figure [5] the dispersion is 
larger for our fiducial model (mi/mj = 0.1, dashed and dot- 
dashed curves) than for a uniform mass function (solid and 
dotted curves). As in Figure |8j the dispersion is smaller for 
the Gaussian profile (dotted and dot-dashed curves) than for 
the uniform source profile (solid and dashed curves) . Perhaps 
the most interesting point is that as the source ellipticity 
increases, the difference between the four curves becomes 
smaller, i.e., the dispersion becomes even less sensitive to 
the mass function and source profile. 
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Figure 11. Dispersion in log versus source size for a positive parity image with different ellipticities and position angles. Rows show 
e=0.3, 0.6, and 0.9, while columns show PA=0, 30, 60 and 90 degrees. Solid and dashed curves have a uniform source profile. The solid 
curve shows a uniform mass function (7711/7712 = 1) and the dashed curve shows a Salpeter mass function with m\/rri2 = 0.1. Dotted and 
dot-dashed curves have a Gaussian profile. The dotted line shows a uniform mass function, while the dot-dashed line shows a Salpeter 
mass function. 



3.5 Accretion Disk Geometry 



Finally, we consider whether variations in accretion disk ge- 
ometry result in observable differences for microlensing. We 
model the source as an annulus with a given half-light radius, 
r$, and hole-to-total area ratio, Q 2 . This approach is useful 
in two ways. First, quasar accretion disks have inner radii 
defined by the innermost stable circular orbit of a particle in 
motion around the central black hole. Second, typical mod- 
els (e.g., the Shakura-Sunyaev disk), emit over a wide range 
of wavelengths. Roughly distinct annular regions within the 
disk are revealed by observations in different bands (see, e.g., 
iMortonson. Schechter fc Wambsgans"sll2005l ). It is important 
to determine whether microlensing can be used to find the 
mass of the central black hole or the scale of the annulus 
within the disk emitting at some wavelength. 



For simplicity we focus on a uniform source. Figure [T21 
shows the dispersion versus source size for disks with Q — 
0.01 (solid curve), 0.1 (dotted), 0.5 (dashed), and 0.9 (dot- 
dashed). For small sources (rs i$ 1.5) the dispersion is nearly 
identical for all values of Q 2 . For larger sources, the dis- 
persion remains similar for Q 2 — 0.01 and for Q 2 = 0.1. 
However, the cases Q 2 = 0.5 and Q 2 = 0.9 have larger dis- 
persions, suggesting that only large holes can significantly 
affect microlensing. 



3.6 Light Curves 

While our analysis has focused on magnification distribu- 
tions, microlensing has a time domain as well and we would 
like to understand whether variability timescales can pro- 
vide more information about the lens and source. Although 
a complete study of microlensing light curves and light curve 
statistics is beyond the scope of this paper, we can exam- 
ine sample light curves and begin to identify useful results. 
We have found that source ellipticity and orientation have 
pronounced effects on the magnification distribution, so we 
now see how they affect light curves. 

To obtain sample light curves, we move the source 
through the magnification map along some trajectory, as 
shown in Figure 1131 The natural time scale is the Einstein 
crossing time, £e = i?rj(m)/«x, where v± is the relative 
transverse velocity between the lens and source. Figure Q3] 
(left) shows the resulting light curves for a source with rs 
— 1 and the same set of ellipticities and orientations used 
in Figure [5] Increasing either the ellipticity or the position 
angle increases the amount of variability, especially on short 
timescales. This is consistent with our previous interpreta- 
tion: small changes to the source position have more effect 
when the source is aligned with the shear (PA=90°) than 
when the source is perpendicular (PA=0°). This can lead 
to a striking amount of rapid variability for highly flattened 
sources. 

To quantify the amount of va ri ability on different 
timescales, we follow iLewis fc Irwiril <|l99rj ) and use the 
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Figure 9. Dispersion in log \fi\ versus source size for a positive 
parity image with different ellipticities and position angles. From 
top to bottom, panels show ellipticities of e=0.3, 0.6, and 0.9. 
The solid, dotted, dashed and dot-dashed curves have position 
angles of PA=0, 30, 60, and 90 degrees, respectively. The source 
is described by a Gaussian profile. The stellar distribution of the 
lens is modeled by a Salpeter mass function with 0.1 < m < 1. 




Figure 10. Illustration of why microlensing magnifications de- 
pend on the orientation of an elliptical source. Both ellipses have 
semimajor axes with lengths a = R^lrh) and ellipticities e = 0.6. 
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structure function as a statistical measure of temporal vari- 
ability. The structure function is defined to be the mean 
square change in the brightness after time At: S(At) = 
{[M{t + At) - M (t)] 2 ) where M is the apparent magnitude 
and the average is over t. To obtain statistically meaningful 
results, we average the structure function over 100 realiza- 
tions of light curves for a given set of parameters. 

The structure functions are shown in Figure [14] (right). 
They all have a roughly linear rise to a plateau beginning 
around At/tE ~ 5. They confirm that there is more vari- 
ability on shorter timescales when the source is elongated 
and/or aligned with the shear. It is customary to define a 
characteristic variability time scale as the interval at which 
the structure funct i on reaches half its p lateau value (see 
iLewis fc Irwm1ll996l ; ISchechter et~aill2003l ) . We see that this 
time scale can vary by a factor of ~2 depending on the el- 
lipticity and orientation of the source. 



Figure 12. Dispersion in log |^| versus source size for a positive 
parity image with different disk geometries. The source is modeled 
as an annulus with a given half-light radius. The solid, dotted, 
dashed and dot-dashed curves have hole-to-total area ratios, Q 2 = 
0.01, 0.1, 0.5 and 0.9, respectively. 



The important implication is that elongated sources 
can lead to more rapid variability than circular sources 
that have identical half-light radii; this effect must be 
taken into account when interpreting observed variability 
timescales. It is not clear wheth er source shape can ex plain 
the rapid variability observed bv lSchechter et all (|2003h and 
IParaficz et al.l (|2006l ) . but it should at least be considered as 
a possible alternative to relativistic motion. 
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Figure 14. Light curves (left) and structure functions (right) for models with e= 0.3, 0.6 and 0.9 (top to bottom). Solid, dotted, dashed 
and dot-dashed curves have PA = 0, 30, 60 and 90 degrees, respectively. 



4 CONCLUSIONS 

We have presented a systematic study of microlensing of fi- 
nite sources. We have extended earlier work by combining an 
finite sour ce with a lens described by a stellar mas s funct ion. 
Following Mor tonson. Schechter fc Wambsganssl (|2005l ) we 
have explored how the surface brightness profile and geome- 
try of the source affect microlensing (JO] and ^3.5|) . We find 
that both effects are of minimal importance although subtle 
differences are apparent: making the source surface bright- 
ness profile steeper (FigHJ) and introducing a large hole in 
the source (Fig |12p both increase the magnification disper- 
sion. 

The mass function can play a more significant role. 
Although the slope of the mass function does not lead 
to noticeable changes in the magnification dispersion, 
the dynamic range can be important. The dispersion for 
a finite source becomes larger as the mass range in- 
creases. This result has been seen before for a broad 
bimodal mass function JSchechter. Wambsganss fc Lewis! 
12004 iLewis fc Gil-Merinoll2006l ), but we have demonstrated 
it for a continuous and relatively narrow mass function as 



would be appropriate for stars. This raises the possibility of 
using microlensing to determine the dynamic range of stellar 
mass functions in distant galaxies. 

Our discussion of dark matter in H3.2I reveals many 
interesting results. We find that the monotonic increase 
in dispersion for a po int source lensed by a mixtur e of 
stars and dark matter (|Schechter fc Wambsganssl l2002h ex- 
tends to the case of combining a small finite source with 
a power-law mass distribution. However, for moderately 
large sources we find that microlensing becomes less pro- 
nounced as the dark mat ter mass fraction is increased. A s 
in previous studies (e.g., ISchechter fc Wambsganssl 120021 ') . 
we find multiple peaks in the magnification histograms for 
moderate dark matter fractions. We also see that nega- 
tive parity images have tails to low magnifications (see 
ISchechter fc Wambsganssl 120021 ') . but only when the source 
is small. 

Finally, we have for the first time considered non- 
circular sources with a range of position angles. We find 
that sources aligned with the shear have larger magnifica- 
tion dispersions than sources orthogonal to the shear. We 
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Figure 13. Construction of light curves. A source with cllipticity 
e = 0.6 and PA= 90° is moved along the trajectory from lower 
left to upper right. The resulting light curves are shown in the 
middle panel of the left column of Figure 14. 



suggest that an aligned source is much more sensitive to the 
position relative to caustics than an orthogonal source. 

Apart from source size, which is fundamentally impor- 
tant, we believe that two of the effects we have identified 
have important physical implications. First, the continuum 
and BLR will be very different in their sensitivity to dark 
matter near a lensed image, particularly a negative-parity 
image. Thus, attempts to measur e the dark matter content 
of ga laxies with microlensing fsee lSchechter fc Wambsganssi 
|2004|) would greatly be nefit from spectroscopic observations 
(see Keeton et alj2006h . Second, elliptical sources, which are 
relevant for inclined disks, may experience much more rapid 
variability than circular sources. This effect will surely be 
important when interpreting microlensing variability time 
scales. 
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